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1 Introduction 



The field of radiative corrections in quantum field theory has always been exciting 
and rapidly developing. At this conference various new methods were presented that 
allow us to keep up with the ever increasing complexity of the problems posed by 
modern particle physics. Some of these methods are concerned with the analj^ic 
evaluation of certain classes of Feynman diagrams (see, e.g., Equally important, 
however, is the development of systematic approximations for complex problems. As 
demonstrated in various physical applications, asymptotic expansions of Feynman 
diagrams prove to be a very efficient tool for this purpose: they provide recipes to 
reduce the number of dimensional scales (masses, momenta) that a diagram depends 
on (see, e.g., 0). The result is a series (possibly asymptotic) in terms of ratios of 
these dimensional parameters. 

However, the validity of an expansion is restricted to a certain - in general finite 
- region of convergence. As we will see. Fade approximations have been used to 
enlarge the validity range of these expansions. In some of the cases described below, 
expansions from different limits could be combined to construct an interpolating 
function which connects the individual, often non-overlapping regions of convergence.]^ 

The outline of this review is as follows: we will begin by describing the general method 
on the basis of the hadronic R ratio. This quantity is extremely important in particle 
physics: not only is it directly measurable at e'^e~ colliders, but it also influences 
other quantities, for example the running of the electro-magnetic coupling constant 
aQEB{s) or the anomalous magnetic moment of the muon. The continuous efforts for 
an accurate evaluation of this quantity have presently reached an accuracy of order 
Og, even though only in the high energy limit. At order a^, due to the successful 
application of the Fade procedure described below, the full energy dependence is 
known. At flrst, the method was applied to non-singlet diagrams which clearly give 
the major contribution to R (cf. Sect.0). Sect.Q will describe the generalizations 
that were necessary to evaluate the singlet contributions. Let us note that recently 
also non-diagonal currents have been taken into account. This opens a new fleld of 
applications related to charged current reactions, like single top production at hadron 
colliders. 

The sections that follow are concerned with a different class of applications of the 
Fade procedure, namely the evaluation of on-shell quantities. First we describe a 
recent calculation of the conversion factor from the quark mass in the MS scheme 
to the on-shell scheme at order a^. Due to the progress in the field of heavy quark 
and top threshold physics, the evaluation of this factor was of utmost importance. 

"'^It might be appropriate to remark that Pade approximations have also been used in the hterature 
to estimate higher order corrections in perturbation theory (e.g. |^). These considerations are of 
completely different nature than the fixed-order predictions which will be discussed in this talk. 
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The second application is the determination of NNLO QCD corrections to the top 
quark decay rate which is a necessary input for a precise experimental determination 
of the top quark properties at future colliders. Closely related is the evaluation of 
the muon decay rate at second order in QED, as well as T{b —>■ ueve) to 0{a1). 
Each of the on-shell quantities above has been calculated by two independent groups 
with complementary methods. The agreement of the results once again confirms the 
validity and accuracy of the Pade method. 



2 General procedure 

We are not going to describe all the details of the procedure for constructing Pade 
approximants, because this has been done in the literature to a sufficient extent (see 
in particular [^,|],^,|^||,^ ) . Nevertheless, for the sake of a closed presentation, let us 
give the main ideas by considering the by now "classic" example of the hadronic R 
ratio: 

, (T(e"''e~ hadrons) . , 

R{s) = ^r^ 1 

cr(e+e /i+yU ) 

The leading two orders in for R{s) are known in analytic form. An analytic 
evaluation of the complete corrections at 0{a1) currently seems to be excluded. The 
exact answer is known only for certain contributions, in particular the terms involving 
a light fermion pair ||10| . 



For the other contributions one has to rely on approximations. There are two obvious 
limiting cases for which this can be achieved: On the one hand, if the center-of-mass 
energy is very large, one may set the quark masses m to zero. One may then use the 
equation 

R{s) = 12n Imn(s/(4m2) + iO+) (2) 

which relates R{s) to the imaginary part of the polarization function 11(2;) along the 
upper branch of the cut z G [1, oo] in the complex plane. Sample diagrams for 11(2;) 
are shown in Fig. |1]. For the moment we will restrict the discussion to diagrams where 
the external currents are connected by a single massive quark line (Fig.|l|(a) and 
(b)). Contributions where the external currents are connected by a massless quark 
line and where the massive quarks couple only to gluons ( "gluon-splitting diagrams" ) 
are numerically unimportant and shall not be addressed here. The modifications for 
diagrams where each of the external currents is connected to a separate Fermion line 
("singlet diagrams", Fig.|I](c)) will be discussed in the next section. 

Taking m = leads to massless propagator diagrams which can be calculated us- 



ing the integration-by-parts algorithm |11| as implemented in the FORM program |12 
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Figure 1: Diagrams contributing to the polarization function n(z). The soHd lines 
are quarks, the springy ones are gluons. The wavy lines represent the external 
currents, (a) and (b) are non-singlet, (c) is a singlet diagram. 
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where the elhpse indicates higher order terms in and in m? /s. 

On the other hand, the leading behavior in the opposite hmit, where the center-of- 
mass energy is close to threshold, i.e. v = a/1 — Am? / s <^ 1, can be deduced from 
the Sommerfeld-Sakharov formula and the two-loop result of the QCD potential (for 
details see 0,0] )• In general, these considerations allow one to deduce the terms that 
are singular for f ^ (cf. Coulomb singularity) as well as the constant term. 

In both limits, however, one can do better. There are well-defined methods to obtain 
expansions around the exact limits m = and s = 4m^ (for reviews see p|JI^). 
These methods reduce the original diagrams that depend on the two scales and 
to single-scale integrals which can be solved analytically. In this way one obtains 
approximations that are valid in certain ranges away from the actual limits. Fig.|^ 
shows the behavior of these expansions as dashed and dotted lines. It immediately 
becomes clear that their validity is restricted to a finite kinematical region. It is the 
purpose of this section to outline the procedure that leads to the solid lines in this 
figure, i.e. the construction of an approximation which is valid over the full v range. 

The limits discussed above are the only two kinematically distinguished points for 
R{s) at 0{a'l) (we disregard the four-particle threshold at s = 16m^ here). Con- 
sidering its connection to the polarization function 11(2;), however, (cf. Eq. (|I|)) there 
clearly is another interesting point, namely z = 0. The coefficients of an expansion of 
Il{z) around z = are called "moments". Because of the cut at 2; = 1, this expansion 
is expected to converge only for \z\ < 1, and at first it seems to be impossible to 
extract any information on R{s) from it. However, one has to recall that 11(2;) is 
analytic everywhere in the complex plane, except along the cut. A mapping of the 
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Figure 2: Threshold expansion up to terms of order vlnv |p!5[ (dotted), high-energy 



expansion up to m /s fTg] (dashed), and Fade result (solid). The singular and 
constant terms in v have been subtracted and the renormalization scale is set to 
fjp = m?. Left: C|-term; Right: CpCA-term (Cp = 4/3 and Ca = 3 are the Casimir 
operators of the fundamental and the adjoint representation of SU(3)). 



form 

- = ^^^ (4) 

transforms the whole z-plane into the unit circle of the tu-plane, such that the upper 
(lower) branch of the cut gets mapped to the upper (lower) semi-circle. The points 
at \z\ = cx) go to = —1, while z = and z = 1 correspond to u = and a; = 1, 
respectively. 

In order to arrive at a smooth function that is easy to approximate, one subtracts 
the threshold singularities from Il{z) (see above) as well as the logarithms of the 
asymptotic high-energy expansion (this has to be done in a suitable way in order not 
to generate logarithms for z ^ 0). After some more manipulations one arrives at a 
function P{uj) whose value at a; = — 1 and whose first few Taylor coefficients around 
a; = are in one-to-one correspondence to the expansion coefficients of n(z) around 
z = and z = oo. 

P{uj) is analytic within |ct;| < 1 and thus its Taylor series around u = converges 
inside this region. However, we need the values of Il{z) on the cut, or equivalently, 
P(|a;| = 1). Thus one has to make sure that convergence of P{uj) is given not only 
for |ti;| < 1, but also for |ti;| = 1. 

This is the point where Fade approximants come into play. The definition of an 
[n/m]-Pade approximant on a function P{uj) is 

ao + aiuj H h anUj"- 
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It has been demonstrated in that such an approximant extends the convergence 
region of the Taylor expansion to its border \uj\ = 1. Performing the mapping back 
to the z-plane, and applying the inverse operations that led from 11(2;) to P{u}), we 
thus have constructed a function that approximates n(z) all over the complex plane, 
including the branches of the cut from z = 1 to z = oo. The imaginary part of this 
function gives rise to the solid lines in Fig.|^. 

A nice feature of Pade approximation is that one has the freedom in varying the 
parameters n and m in Eq. (^). If convergence of the Pade approximants was not 
given, this would manifest itself in strong variations of the result for different values 
in n and m. For 11(2;), for example, this dependence is so weak that different Pade 
approximants produce curves that would be hardly distinguishable from one another 
in Fig.|(see, e.g., i). 

As it was mentioned above, this procedure was applied for the first time in [^] to 
the three-loop QED vacuum polarization. It was later on generalized to the QCD 
case [|] for various external currents (vector, axial- vector, scalar, pseudo-scalar). 
In all of these papers, only the leading two terms in the asymptotic expansion around 
z oo were taken into account. After higher order terms in this limit became 
available [0, the method was extended to include them |^ and a further stabilization 
of the Pade predictions was observed [p!7||9|. Recently the method has been apphed 
to non-diagonal currents in order to derive the dominant contributions to single top 
production at hadron colliders. 



3 Pade approximation for singlet diagrams 

Above we described in some sense the "optimal case" : information on both the limits 
z ^ and z ^ oo was available, and the leading threshold behavior was known. 
Furthermore, the analytic structure of 11(2;) was such that the expansion around z = 
had the form of a plain Taylor expansion as opposed to an asymptotic series, i.e., it 
did not contain any logarithms of z. This is due to the fact that the discussion was 
restricted to the non-singlet contributions. The corresponding diagrams (cf. Fig.|^(a) 
and (b)) do not have massless cuts: the external currents are connected by a single 
massive quark line, and cutting the diagram in halves always involves a cut through 
this line. 

In the remaining part of this review we will be concerned with exceptions to this 
"optimal procedure." The first case we will consider are the singlet contributions 
to 11(2;). They are distinguished from the non-singlet diagrams in the sense that 
the external currents are each connected to separate Fermion lines which in turn are 
connected to each other by gluons (cf. Fig. ||(c)). If the external currents are vector- 
like, these diagrams vanish at 0{a1) due to Furry's theorem. If, however, one is 
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concerned with axial-vector, scalar, or pseudo-scalar currents, massless cuts occur.0 

These massless cuts spoil the analyticity of P{uj) (see Sect.^ within < 1, and 
thus also the convergence of its expansion around u = 0. Luckily, for the singlet 
diagrams the analytic expressions of the massless cuts are known. Thus, denoting 
these massless cuts as -Rmi(s), we may employ the dispersion relation to write 



Us{z) =U^i{z) +fl{z) , with Il{z)=C-^ ds ^ 



oo 



and n^i(z) = C-^ f ds ^^^^\ = C-^ f ds 
Jo s- Am^z Jo 



' (6) 



s — Airfiz 



{C = 1271 for external vector and axial-vector currents, C = Sn for scalar and pseudo- 
scalar currents). Here, ns(-2;) is the singlet contribution to the polarization function 
and 

Rs{s) = ClmUs{s/{4m^) + zO+) . (7) 

tl{z) is analytic in the complex plane cut along z G [l,oo]. It can be obtained by 
evaluating Il^\{z) through Eq. and subtracting it from lis (2;), i-e., the result for 
the singlet diagrams (see Fig.[l| (c)). One can then apply the Fade procedure outlined 
in Sect.|^ to Il{z). 

For details on the evaluation of the integral for nmi(2;) in Eq. (^) and the results we 
refer to [§]. At this point it shall be sufficient to mention that the combination of 0,|19 



(non-singlet) and (singlet) provides the current knowledge of R{s) at C(a^). Some 
contributions are known analytically, and the accuracy of the approximations in all 
the other cases is extremely good. Let us also remark that at 0{al), Il{z) is known 
in its high-energy expansion up to the terms oc m'^/s'^ |20, 2l] , p2| . No moments for 



0—^0 are available yet, and therefore a Fade approximation along the lines of the 
previous section is still out of reach. 

As another application of the methods described above let us note that the results at 
(9(a^) were combined with the one-loop electro- weak corrections in order to predict 
the total cross section for e~^e~ — tt at a linear collider [ESl. 



4 Relation between MS and on-shell quark mass 

The MS scheme is a very convenient renormalization scheme, in particular from the 
technical point of view. Renormalization constants in the MS scheme do not depend 

^ In the axial- vector case, the purely gluonic cuts are zero according to the Landau- Yang theorem; 
but in order to avoid the axial anomaly, it is convenient to consider a full SU(2) doublet like {t, b). 
Taking mi, = 0, the massless cuts arise from cuts involving b quarks. 
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Figure 3: Diagrams contributing to (a) Zm, (b) T{t bW), and (c) r(/x — > euei^fj,)- 



on any masses or momenta, which means that their evaluation provides a certain 
freedom in choosing the particular set of diagrams to calculate. For example, the 
four-loop quark anomalous dimension was computed using two completely different 
approaches . 



However, comparison of the theoretical results to experimental data often requires 
to express the involved quantities in terms of their on-shell values. This is why the 
conversion factor 

Zm = — {^) 

m 

that relates the on-shell to the MS quark mass is of great importance. In order to 
obtain this quantity, one has to evaluate the quark propagator 

E(g)=^Sv(g')+mSs(g^) (9) 

at the on-shell point = m^, where q is the external momentum and m is the quark 
mass. A sample diagram that contributes to Tj{q) at order is shown in Fig.|^(a). 

Technically, these diagrams carry only a single scale and should be accessible through 



the integration-by-parts algorithm |TT[. However, it turns out that the level of com- 
plexity for n-loop on-shell diagrams is comparable to (n + l)-loop massive tadpole 
diagrams, for example. At the two-loop level, the problem of finding the recurrence 



relations was solved in 25 



The three-loop case seemed to be out of reach for quite some time. This is why the 
first calculation of the conversion factor was performed using a different, semi-analytic 



approach with the help of Fade approximants |^6[. Looking at the Fermion propagator 
as a function of 2; = g^/m^ in the complex z plane, we find a similar structure as for 
the polarization function Il{z) in Sect.|^. Es v(g^) is analytic in the complex plane 
cut along z & [1, 00]. Thus, in principle one can follow the same strategy as for H(z): 
based on the expansions around = and q^ = oc one constructs an approximation 
for Es,v('?^) in the whole q^ plane, including q^ = m? (2; = 1), the point of interest. 

A complication one has to face here is that the functions T,^y{q^) depend on the 
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Figure 4: Gauge parameter dependence of —g{z). {g{z) is constructed from TiQ,y{q^) 
in such a way that g{l) = —Zm — for details see |2^). "F" denotes the 0{as) 
contribution, "FF" the C| term at 0{a1). The sohd hne corresponds to ^ = —2, the 
others to ^ = —5, 0, +2, +5 (from short to long dashes) [Figures taken from p^j. 



strong gauge parameter ^ in general.^ Only at the on-shell point = m? are they 
gauge independent. Thus, the expansions around 2; = and 2; = 00 exphcitly contain 
^. If they could be re-summed exactly, ^ would drop out in the limit z — > 1. But here 
the exact re-summation will be replaced by a Fade approximation and thus the result 
depends on even for z = 1. However, the argumentation is that the dependence 
on ^ is weak at z = 1, in the sense that any "reasonable" choice of ^ leads to valid 
predictions, with the spread among different choices being within the error of the 
Fade approximation. The question of what a "reasonable" choice is can be answered 
by making the natural assumption that a smooth curve gets better reproduced by 
Fade approximants than a strongly varying one. Looking at Fig.§, it is not hard to 
decide that the favored values for ^ are within a few units around ^ = 0. 

Following the outlined procedure, the authors of were able to deduce the value 
of Zm with an uncertainty of around ±3%. Shortly after this result was presented, a 
second group published the analytical result for the conversion factor Zm [^. They 



managed to establish the recurrence relations for three-loop on-shell diagrams derived 
by integration-by-parts, which provides an extremely useful tool for various other 
applications (see also |2^). 



^We define the gluon propagator in i?^ gauge as i{—g^i, + ^q^iQu / q^) / {q^ + «0+). 
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5 Top decay to and second order QED corrections to muon 
decay 

The decay rate of the top quark is expected to be measured at future coUiders with 
about 10% accuracy. This is roughly the size of the C(a<j) corrections to this quantity, 
which is why the evaluation of the contribution was necessary. One option to 
evaluate this decay rate is to calculate the on-shell top quark propagator to the 
appropriate order and to take the imaginary part. A typical diagram whose imaginary 
part contributes to the rate r{t bW) at order is shown in Fig.|^(b). In a first 
approximation, one may take the bottom quark and the W boson to be massless. Then 
the lowest order diagram is a massless propagator. Starting from 0{as), however, 
one is again faced with proper on-shell diagrams. In contrast to the problem of the 
MS/on-shell conversion factor (see Sect.^ the diagrams for t —>■ bW contain massless 
cuts due to the presence of the b and the W. The solution of p9[ (see also |Q) to this 
problem was to evaluate the expansion around /m^ = and to take the imaginary 
part before evaluating a Fade approximation. It is therefore not possible to take 
information from the limit q"^ —>■ oo into account, because otherwise one receives 
contributions to the imaginary part coming from the cut starting at = mf. 

Nevertheless, the Fade approximants constructed from the expansion around = 
alone give a fairly accurate result (judging from the spread of the different approx- 
imants). In addition, it agrees well with an earlier result which relied on an 
expansion around 1 — ml/mf = 0. Considering the fact that these approaches are 
based on expansions around two completely different limits - which, in addition, are 
both far from the physical point - their agreement within very small error bars is a 
clear demonstration of the power of the applied methods. 

Another example of this kind is the decay rate of the muon. Following a strategy 
closely related to the one described above for top decay, it was possible to evaluate 
the second order QED corrections to this quantity in a semi-analytical way ||32|| . 
Note, however, that the corresponding diagrams contain four closed loops here (cf. 
Fig.^(c)), even though one of them (made up by the neutrino lines) is always a 
massless self-energy insertion. Also in this case the Fade approximants agree nicely 



with the previously known analytical result ||33[ which provides an important check 
on the latter. 

Along the same lines one can also obtain the corrections of order a"^ to the decay rate 
r{b — s> ueue), the only technical difference being the presence of non-Abelian diagrams. 
Also here the semi-analytical result |^ approximates the analytical formula to 
high accuracy. 
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6 Conclusions 



We reviewed the method and recent apphcations of Pade approximation to fixed order 
calculations in QCD. Originally developed for the hadronic e+e~ cross section, the 
approach proved to be useful also for a completely different class of problems related 
to on-shell phenomena, for example in heavy quark physics. Let us conclude by 
pointing out that the continuously refining field of expansion techniques for Feynman 
diagrams (see [0]) should pave the way to numerous new applications for the Pade 
method. 
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